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Stau relic density at the Big-Bang nucleosynthesis era consistent with the abundance 
of the light element nuclei in the coannihilation scenario 
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We calculate the relic density of stau at the beginning of the Big-Bang Nucleosynthesis (BBN) era 
in the coannihilation scenario of minimal supersymmetric standard model (MSSM). In this scenario, 
stau can be long-lived and form bound states with nuclei. We put constraints on the parameter 
space of MSSM by connecting the calculation of the relic density of stau to the observation of the 
light elements abundance, which strongly depends on the relic density of stau. Consistency between 
the theoretical prediction and the observational result, both of the dark matter abundance and the 
light elements abundance, requires the mass difference between the lighter stau and the lightest 
neutralino to be around lOOMeV, the stau mass to be 300 - 400 GeV, and the mixing angle of the 
left and right-handed staus to be sin# T = (0.65 - 1). 
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I. INTRODUCTION 

Cosmological observations have established the exis- 
tence of the non-baryonic dark matter (DM) [l| . These 
observations suggest that the DM is a stable and weakly- 
interacting particle with a mass of 0(100) GeV. Many 
hypothetical candidates for the DM have been proposed 
in models of particle physics beyond the standard model 
(SM), and one of the most attractive candidates is the 
lightest neutralino, x°, in supersymmetric extensions of 
the SM with R parity conservation. Neutralino is a linear 
combination of the superpartners of U(l), SU(2) gauge 
bosons and the two neutral Higgses, and is stable when 
it is the lightest supersymmetric particle (LSP). Indeed 
it accounts for the observed DM abundance when it is 
degenerate in mass to the next lightest supersymmetric 
particle (NLSP) and hence coannihilates with the NLSP 
We consider the setup that the LSP is a neutralino 
consisting of mainly bino, the superpartner of U(l) gauge 
boson, and the NLSP is the lighter stau, the superpart- 
ner of tau lepton. This is naturally realized in the MSSM 
with the unification condition at the grand unified the- 
ory scale. The minimal supersymmetric SM (MSSM) has 
two eigenstates of stau as physical state. In absence of 
inter-generational mixing the mass eigenstate of stau is 
given by the linear combination of the left-handed stau 
fh and the right-handed stau tr as 
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where 9 T is the left-right mixing angle and 7 T is the CP 
violating phase. 

In a scenario of the coannihilation, the NLSP stau can 
be long-lived if the mass difference, 8m, between neu- 
tralino and stau is small enough to forbid two-body de- 
cays of stau into neutralino. It was shown in Q that the 



lifetime of stau is longer than 1000 second for 8m < 100 
MeV. It is known [4H18I] that the long-lived charged par- 
ticles affect the relic abundance of the light nuclei during 
or after the big-bang nucleosynthesis (BBN). 

There is a discrepancy, the so-called 7 Li problem, on 
the primordial 7 Li abundance between the prediction 
from the standard BBN (SBBN) and the observations 
[]~9l . |20| . Combined with the up-to-date values of baryon- 
to-photon ratio, v = (6.225±0.170) x 10~ 10 from Wilkin- 
son Microwave Anisotropy Probe (WMAP) [1[, the 
SBBN predicts the 7 Li to proton ratio, ( 7 Li/H)sBBN = 
5.24^Q gy x 10~ 10 , which is by about four-times larger 
than its observed value in poor- metal halos [19|, l2l| . Be- 
cause there exists no general agreements about astro- 
physical scenarios to reduce the 7 Li abundance [22T4251 ] . 
it is natural to consider nonstandard effects. 

The authors have investigated the BBN including the 
long-lived stau [3, EH- The long-lived stau form a 
bound state with nuclei {fN), and consequently con- 
vert it into a nucleus with a smaller atomic number. 
Here N stands for a nucleus. In this scenario, the abun- 
dance of 7 Li is reduced through the conversion process, 
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- 7 Li and the further destruction of 7 Li 



by either a collision with a background proton or another 
conversion process, (f 7 Li) — > x° + v t + 7 He. Therefore, 
the more these bound states are formed, the more the 
7 Li abundance is reduced. The number density of the 
bound state is determined by the relic density of stau. In 
[01 [HI j we assumed Yf.po and 8m to be free parameters 
of the scenario, where Y? is the yield value of stau 
at the time of decoupling from the thermal bath. The 
full calculation of the light nucleus abundances showed 
a region in {8m, Yf^o) plane where the 7 Li problem is 
solved consistently with the observational constraints on 
the other nuclei. The region points If fo to be close to 
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the yeild value of the DM and 8m to be around 0.1 GeV, 
respectively. 

In this work, we improve our previous analysis by cal- 
culating I^fo with taking the relic abundance of DM 
into account. The outline of this paper is as follows. In 
section [ill we review the formalism for the calculation 
of the relic density of stau, and derive the Boltzmann 
equations of stau and neutralino for the calculation. In 
section IIII1 we present the numerical results of the cal- 
culation, and prove the parameter space for solving the 
7 Li problem. Section IIVI is devoted to a summary and 
discussion. 



II. FORMALISM FOR THE CALCULATION OF 
RELIC DENSITY OF STAU AT THE BBN ERA 

In this section, we prepare to calculate the relic density 
of stau at the BBN era. Firstly, in subsection III Al we 
briefly review the Boltzmann equations for the number 
density of stau and neutralino based on the thermal relic 
scenario. Then, in subsection lll B[ we discuss the number 
density evolution of stau and neutralino quantitatively. 
In subsection lll CI we investigate the significant processes 
for the calculation of the relic density of stau. Finally, 
we obtain the Boltzmann equations for the relic density 
of stau in a convenient form. 



Boltzmann equations for the number density 
evolution of stau and neutralino 



In this subsection, we show the Boltzmann equations of 
stau and neutralino and briefly review their quantitative 
structure based on [2^| (see also a recent paper 27]). 

We are interested in the relic density of stau in the 
coannihilation scenario. In this scenario, stau and neu- 
tralino are quasi-degenerate in mass and decouple from 
the thermal bath almost at the same time @. Thus the 
relic density of stau is given by solving a coupled set of 
the Boltzmann equations for stau and neutralino as si- 
multaneous differential equation. For simplicity, we use 
the Maxwell-Boltzmann statistics for all species instead 
of the Fermi-Dirac for fermions and the Bose-Einstein 
for bosons, and assume T invariance. With these simpli- 
fications, the Boltzmann equations of them are given as 
follows 
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Here n and n eq represent the actual number density and 
the equilibrium number density of each particle, and H 
is the Hubble expansion rate. Index i denotes stau and 
neutralino, and indices X and Y denote SM particles. 
Note that if relevant SM particles are in thermal equi- 
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then these equations are reduced into a familiar form. 
(av) and (a'v) are the thermal averaged cross sections, 
which is defined by 
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where / is the distribution function of a particle, v is 
the relative velocity between initial state particles, and 
g 12 = 2(1) for same (different) particles 1 and 2. In this 
work, we assume that all of the supersymmetric particles 
except for stau and neutralino are heavy, and therefore 
do not involve them in the coannihilation processes. 

The first line on the right-hand side of Eqs. ((5J), ([3]), 
and ([4J accounts for the annihilation and the inverse 
annihilation processes of the supersymmetric particles 
(ij -H> XY). Here index j denotes stau and neutralino. 
As long as the R-parity is conserved, as shown later, 
the final number density of neutralino DM is controlled 
only by these processes. The second line accounts for 
the exchange processes by scattering off the cosmic ther- 
mal background (iX O jY). These processes exchange 
stau with neutralino and vice versa, and thermalize them. 
Consequently, the number density ratio between them is 
controlled by these processes. Instead, these processes 
leave the total number density of the supersymmetric 
I particles. Note that in general, although there are terms 
>which account for decay and inverse decay processes of 
I stau (f O xXY...) in the Boltzmann equations, we omit 
them. It is because we are interested in solving the 7 Li 
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problem by the long-lived stau, and the whole intention 
of this work is to search parameters which can provide 
the solution for the 7 Li problem. Hence we assume that 
the stau is stable enough to survire until the BBN era, 
and focusing on the mass difference between stau and 
ncutralino is small enough to make it possible. 

B. The evolution of the number density of stau 
and neutralino 

In this subsection, we discuss the evolution of the num- 
ber density of each species. Firstly, we discuss the num- 
ber density evolution of neutralino DM. Since we have as- 
sumed R-parity conservation, all of the supersymmetric 
particles eventually decay into the LSP neutralino. Thus 
its final number density is simply described by the sum 
of the number density of all the supersymmetric particles 



This is because the cross sections of the exchange pro- 
cesses are in the same order of magnitude as that of the 
annihilation and the inverse annihilation, but the num- 
ber density of the SM particles is much larger than that 
of the supersymmetric particles which is suppressed by 
the Boltzmann factor. Thus even if the total number 
density of stau and neutralino is frozen out at the tem- 
perature Tf, each number density of them continue to 
evolve through the exchange processes. 

Thus, to calculate the relic density of stau, we have 
to follow the two-step procedures. As a first step, we 
calculate the total relic density of the supersymmetric 
particles by solving the Eq. (J7J). We use the publicly 
available program micrOMEGAs |28[ to calculate it. The 
second step is the calculation of the number density ratio 
of stau and neutralino. The second step is significant for 
calculating the relic density of stau at the BBN era, and 
hence we will discuss it in detail in the next subsection. 



(6) 



For N, that is the number density of the neutralino, we 
get the Boltzmann equation by summing up Eqs. ([2]). 
O, and Q, 
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Notice that the terms describing the exchange processes 
in each Boltzmann equations cancel each other out. Solv- 
ing the Eq. (J7J, we obtain N and find the freeze out 
temperature of the total number density of all the super- 
symmetric particles Tf by using the standard technique 
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Here, g and are the internal degrees of freedom and 
the mass of neutralino, respectively. The Planck mass 
iripi = 1.22 x I0 19 GeV, and g* are the total number of 
the relativistic degrees of freedom. Consequently, we see 
that 4 GeV < T f < 40 GeV for 100 GeV < m s < 1000 
GeV. 

Now we will discuss the number density evolution of 
stau. To obtain the relic density of stau, we solve a cou- 
pled set of the Boltzmann equations, ([2]), (j3j, and (j4j, 
as simultaneous differential equation. Each Boltzmann 
equation contains the contributions of the exchange pro- 
cesses. These processes exchange stau with neutralino 
and vice versa. At the temperature Tf, the interaction 
rate of the exchange processes is much larger than that of 
the annihilation and the inverse annihilation processes. 



C. The exchange processes and Lagrangian for 
describing them 

After the freeze-out of the total number density of stau 
and neutralino, each of them is exchanged through the 
following processes 
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Notice that although there are other exchange processes 
via weak interaction (for example, tW -f-> x^r, tu t *^ 
)(W, and so on), we can omit them. This is because 
the number density of W boson is not enough to work 
these processes sufficiently due to the Boltzmann factor 
suppression, and the final state W boson is kinematically 
forbidden when the thermal bath temperature is less than 
Tf. These processes (Eq. (TlU)) ) are described by the 
Lagrangian 
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(11) 



where e is the electromagnetic coupling constant, Pl and 
Pr are the projection operators, and I € {e,fj,}. gh and 
gn are the coupling constants, given by 
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where g is the SU(2)l gauge coupling constant, and 9w 
is the Weinberg angle. 

The evolution of the stau number density is governed 
only by the exchange processes fEq.(fl"0|)') after the freeze- 
out of the total relic density of stau and neutralino. 
When we calculate it, we should pay attention to two 
essential points relevant to the exchange processes. 
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FIG. 1: The inverse annihilation rate of tau leptons (r)n^, 
the inverse decay rate of tau leptons (T), and the Hubble 
expansion rate H as a function of the thermal bath tempera- 
ture. 



of tau lepton. When the SM particles X, Y, and 



Z are in the thermal equilibrium, (nyriz) / '(n^n^) = 
[nxny ...) I (ji^ny ...) = 1, and hence tau leptons are 
sufficiently produced through the inverse annihilation 
and/or the inverse decay processes as long as these inter- 
action rates are larger than the Hubble expansion rate. 
Therefore, whether tau leptons are in the thermal bath 
or not can be distinguished by comparing the Hubble ex- 
pansion rate H with the inverse annihilation rate of tau 
lepton (av)n e ^, and the inverse decay rate of tau lepton 
(r). In other words, the inequality expression 



(av)n eq 



x 



> H and/or (T) > H (16) 



indicates that tau leptons are in the thermal bath. Fig. 
1 shows (r)n^, (r), and H as a function of the thermal 
bath temperature. As shown in Fig. 1, the inverse decay 
rate of tau lepton is much larger than the Hubble expan- 
sion rate. Thus, we can conclude that tau leptons remain 
in the thermal bath still at the beginning of the BBN. 



One is the competition between the interaction rate of 
the exchange processes and the Hubble expansion rate, 
since when these interaction rates get smaller than the 
Hubble expansion rate, the relic density of stau would be 
frozen out. The other is whether tau leptons are in the 
thermal bath or not. The interaction rate of the exchange 
processes strongly depends on the number density of tau 
leptons. When tau leptons are in the thermal bath, the 
number density ratio between stau and neutralino are 
given by the thermal ratio, 
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D. Calculation of the number density ratio of stau 
and neutralino 



We are now in a position to calculate the number den- 
sity ratio of stau and neutralino. In this subsection, we 
will show a set of relevant Boltzmann equations. 

The right-hand side of the Boltzmann equations (Eqs. 
©, (O, and flU) depends only on temperature, and 
hence it is convenient to use temperature T instead of 
time t as independent variable. To do this, we reformu- 
late the Boltzmann equations by using the ratio of the 
number density to the entropy density s : 



through the exchange processes. On the contrary, once 
tau leptons decouple from the thermal bath, the ratio 
cannot reach this value. To calculate the relic density 
of stau, we have to comprehend the temperature of tau 
lepton decoupling. 

To see whether tau leptons are in the thermal bath or 
not, we consider the Boltzmann equation for its number 
density, n T , 
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Consequently, we obtain the Boltzmann equations for the 
number density evolution of stau and neutralino 
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where indices X, Y, and Z denote the SM particles, 
and (r) represents the thermal averaged decay rate 
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Here 3* (T) is the relativistic degrees of freedom, and we 
use 

s = ^g*(T)T 3 , H = IMgV 2 — . (21) 
45 rripi 

We obtain the relic density of stau at the BBN era by 
integrating these equations from Tf to the temperature 
for beginning the BBN under the initial condition of the 
total number density of stau and neutralino. These equa- 
tions make it clear that if the tau number density is out 
of the equilibrium, the ratio between those of stau and 
neutralino does not satisfy the Eq.(TT3|). 



III. NUMERICAL RESULTS 
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FIG. 2: Total abundance of staus and neutralinos, which cor- 
responds to the relic abundance of DM. Each line shows the 
total abundance, and attached value represents the stau mass. 
Yellow band represents the allowed region from the WMAP 
observation at the 3<r level , and the region inside the horizon- 
tal dotted lines corresponds to allowed region at the 2a level 
[J. In the left side of vertical lines, the LSP is left-handed 
sneutrino. Three lines correspond to m? = 400, 350, 300 GeV 
from left to right. 



In this section, we will first show the evolution of the 
stau number density, and then study the relation between 
the relic density of stau and the modification of nucle- 
osynthesis. Finally, we study a solution of the 7 Li prob- 
lem with long-lived stau in the coannihilation scenario 
based on Ref. ^M^- 



A. Total abundance 

As a first step for the calculation of the relic number 
density of stau, based on the discussion in section lllBl we 
calculate the total abundance of stau and neutralino with 
micrOMEGAs [28]. Fig. [5] shows the total abundance, 
which corresponds to the relic abundance of DM, as a 
function of sin0 r , where 9 T is the mixing angle between 
left and right-handed stau. Each curved line shows the 
total abundance for each stau mass, and horizontal band 
represents the allowed region from the WMAP observa- 
tion at the 3ct level (0.0913 < fl DM h 2 < 0.1285), and 
the region inside the horizontal dotted lines corresponds 
to allowed region at the 2a level (0.0963 < flumh 2 < 
0.1213) 11 In the left side of vertical lines, the LSP 
is the left-handed sneutrino. Three lines correspond to 
m? = 400, 350, 300 GeV from left to right. Since the 
left-handed sneutrino DM has been ruled out by con- 
straints from the direct detection experiments [29j, we 
focus on the right-side region. Here we took 7 T = and 
Sm = 100 MeV. 

The total abundance increases first as the heavier stau 
mixes to the lighter stau, then turns to decrease at 



sia0 T ~ 0.8. The increase of the abundance can be un- 
derstood by the fact that the annihilation cross section of 
f + f — > r + r becomes smaller as the heavier stau mixes. 
The increase is gradually compensated by two annihila- 
tion processes, f+f* -> W + +W~ andxi+T — > W~+v T , 
as the left-right mixing becomes large. The latter process 
can not be ignored because left-handed sneutrino is de- 
generate to stau and neutralino in the present parameter 
set. These processes become significant for sin^,- < 0.8. 
Another process which reduces the total abundance due 
to the left-right mixing is f+f* —> t+i through s-channel 
exchange of the heavy Higgses. This annihilation process 
becomes significant as the mixing reaches to 7r/4 and the 
masses of two staus are split. In fig. 2, the mass differ- 
ence between the lighter and the heavier stau is fixed to 
be 30 GeV to maximize the DM abundance, but the same 
result can be obtained by changing the stau mass for an- 
other values of the mass difference. As shown in Fig. [21 
the total abundance also strongly depends on m?. This is 
understood as follows. In the non-relativistic limit, since 
the relic number density of relic species is proportional to 
(m re iic{av) S nr S )~ 1 and (av) sum (Eq. (JgJ) is proportional 
to 1/m-re/ic I2fl|> the total number density N is propor- 
tional to rrif, 
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and the total abundance is given by QnMh 2 ~ m?N. 
Thus the total abundance is proportinal to to?, and it is 
consistent with the result in Fig. [5J 
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FIG. 3: The evolution of the number density of negative 
charged stau. Each line attached [8m] shows the actual evo- 
lution of the number density of stau, while the one atattched 
[8m (thermal)] shows its evolution under the equilibrium de- 
termined given by Eq. ()13|) and the total relic abundance. 
Yellow band represents the allowed region from the WMAP 
observation at the 2a level 



B. Stau relic density at the BBN era 



Next, we solve the Boltzmann equations (|T8|) . (|19l) . 
and (|20l) numerically, and obtain the ratio of the stau 
number density to the total number density of stau and 
neutralino. Fig. [3J shows the evolution of the number 
density of stau as a function of the universe tempera- 
ture. Here we took m? = 350 GeV, sin6> T = 0.8, and 7 T 
= and chose 8m = 10 MeV, 50 MeV, and 100 MeV as 
sample points. Each line attached [8m] shows the actual 
evolution of the number density of stau, while the one 
atattched [Jm(thermal)] shows its evolution under the 
equilibrium determined by Eq. (I13|) and the total relic 
abundance. Horizontal dotted line represents the relic 
density of DM, which is the total abundance calculated 
above. We took it as a initial condition of total value for 
the calulation of the number density ratio. Yellow band 
represents the allowed region from the WMAP observa- 
tion at the 2cr level [![. 

The number density evolution of stau is qualitatively 
understood as follows. As shown in Fig. [31 the freeze- 
out temperature of stau almost does not depend on 8m. 
It is determined by the exchange processes Eq. (|T0| . 
whose magnitude (av}Y?Y~ is governed by the factor 
e -(m T -8m)/T ^ wuere TOt represents the tau lepton mass. 
The freeze-out temperature of the stau density T/( ratio ) is 
given by (m T — <5m)/T/( ra tj ) — 25 as in Eq. ©, since the 
cross section of the exchange process is of the same mag- 
nitude as weak processes. Thus T/( ra tio) hardly depends 
on 8m. In contrast, the ratio of the number density be- 
tween stau and neutralino depends on 8m according to 



Eq. (IT3l) . rif/n-x ~ exp(— 8m/T), since they follow the 
Boltzmann distribution before their freeze-out. Thus, the 
relic density of stau strongly depends on 8m. 

Here, we comment on the dependence of the stau relic 
density rif- on other parameters such as mf, 6 T , and 
7 T . The number density of the negatively charged stau 
is expressed in terms of the total relic density N by 



N 



2(1 



Jm/Tf 
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Here, the freeze-out temperature J/( rat j ) hardly depends 
on these pararameters. This is because the cross sec- 
tion of the exchange processes are changed by these pa- 
rameters at most by factors but not by orders of magni- 
tudes, and the Tj( rat i D ) depends logarithmically on (an) 
as shown in Eq. ([5]). On the other hand, the total relic 
density N is proportional to mf as in Eq. ([22]) . The 
value of N is also affected by the left-right mixing 6 T 
as seen in Fig. [2] since the annihilation cross section de- 
pends on this parameter. In contrast, j T scarcely affects 
the relic density, since this parameter appears in the an- 
nihilation section through the cross terms of the contri- 
butions from the left-handed stau and the right-handed 
one, and such terms always accompany the suppression 
factor of m T /mf compared to the leading contribution. 
Thus the relic number density of stau rif strongly de- 
pends on mf and 9 T while scarecely depends on j T . 

We comment on the generality of our method to calcu- 
late the density of exotic heavy particles that coannihi- 
late with other (quasi)stable particles: we calculate the 
total number density of these particles and then calcu- 
late the ratio among them by evaluating the exchange 
processes such as Eq. (fTOf . This method of calculation 
can be found versatile in various scenarios including the 
catalyzed BBN and the exotic cosmological structure for- 
mation 13114331 . 



C. Long-lived stau and BBN 

After the number density of stau freezes out, stau de- 
cays according to its lifetime 3], or forms a bound state 
with a nuclei in the BBN era. Their formation rate has 
been studied in literatures [3, EH, E3] • The bound states 
modify the predictions of SBBN, and make it possible to 
solve the 7 Li problem via internal conversion processes 
in the bound state [l4 - |l6f . 

In Fig. [4] we show parameter regions that are consis- 
tent with the observed abundances of the DM and of the 
light elements. We calculate the relic density of stau by 
varying the value of 8m with the values of mf = 350 
GeV, sin^T = 0.8, and 7 T = 0. With these parameters, 
the allowed region is shown inside the dotted oval. We 
see that there are allowed regions at Y f - ~ 10~ 13 - 10~ 12 
for 8m < 130 MeV to solve the 7 Li problem at 3d. On 
the other hand, it is found that the observational 6 Li to 
7 Li ratio excludes Y f - > 10" 15 and 8m < 100 MeV. We 
will explain this feature as follows. 
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FIG. 4: Parameter regions that are consistent with the ob- 
served abundances of the DM and the light elements. Con- 
straints from the observed 7 Li abundances are due to Bonifa- 
cio et al. (3(J and Melendez and Ramilez [U ■ The calculated 
relic number density of stau with the indicated parameter is 
also shown. The top left region is excluded from 6 Li overpro- 
duction for K- > 1(T 15 and 8m < 100 MeV. 



We have adopted following observatinal abundances of 
6 Li and 7 Li. Throughout this subsection, ni denotes the 
number density of a particle "z" , and observational errors 
are given at la. For the ne,^ to m L i ratio, we use the 
upper bound [34j . 



(n«Li/«'Li)„ < 0.046 ± 0.022 + 0.106, 



(24) 



with a conservative systematic error (+0.106) [Hj]. For 
the Li abundance, we adopt two observational values of 
the n.7Li to n B ratio. Recently it has been reported to be 



logio(n7 U /n H )j ; 



-9.90 ±0.09, 



(25) 



by Ref. [361 ] . and on the other hand, a milder one was 
also given by Ref. [37J , 



log 10 (m Li /n H ) P = -9.63 ± 0.06. (26) 

In the current scenario 6 Li can be overproduced by the 
scattering of the bound state ( 4 Hef ~) off the background 
deuterium through ( 4 Hef~) + D — > 6 Li + f~ 4]. The 
abundance of the nonthermally-produced 6 Li through 
this process is approximately represented by 



with (c7"y)eLi the thermal average of the cross section 
times the relative velocity for this process [Hj], and n B 
the number density of deuterium. By using ([24]) . we 
see that the additional 6 Li production is constrained to 
be AYe L i < 0(1O -21 ). Numerical value of {av)^ gives 
(av)6 U n D /H ~ 0(1O- 6 ) at T -10 keV. Then from (27|) 
it is found that the upper bound on the abundance of 
stau should be Yf- < 10~ 15 . Because the bound state 
( 4 Hef~) forms at T <10 keV, this process is strongly 
constrained for t>- > 10 4 with r f - being the stau life- 
time, which corresponds to 6m < 100 MeV. Note that 
the ratio (av)ei J - 1 nD / H rapidly decreases as the cosmic 
temperature decreases, and this nonthermal production 
of 6 Li is much more effective just after formation of the 
bound state. This is the reason why we can estimate ([27]) 
at around 10 keV. 

On the other hand, the rates of 7 Be and 7 Li destruc- 
tion through the internal conversion [l4T - tl6l ] could be 
nearly equal to the formation rates of the bound sate 
( 7 Bef~) and ( 7 Lif~), respectively. This is because the 
timescale of the destruction through the internal conver- 
sion is much faster than that of any other nuclear reaction 
rates and the Hubble expansion rate. Then the destroyed 
amount of 7 Be (or 7 Li after its electron capture) is ap- 
proximately represented by 



(g-f)bnd,7^ Bc 

H 



(28) 



(27) 



where (av) hnd ,7 ~ lO^GeV^CT/SOkeV)- 1 / 2 ^) 2 x 
(A/7)- 3 / 2 (£ 6 7 Bc /1350keV) is the thermally -averaged 
cross section times the relative velocity of the bound- 
state formation for ( 7 Bef~) 0, We request AYr Bc 
to become — 0(1O~ 20 ) to reduce the abundance of 7 Be 
to fit the observational data. Then the abundance of 
f ~ should be the order of AV7 Be ((cn;)bnd,7W Be / ' H)^ 1 — 
O(10~ 12 ) with Mbnd,7m Bo /tf ~ 10" 8 at T = 30 keV. 
Because (o"i>)bnd,7 m Bc/^ decreases as the cosmic tem- 
perature decreases (a T 1 / 2 ), the destruction is more ef- 
fective just after the formation of The bound state. This 
validates that we have estimated (|2"51) at 30 keV. There- 
fore the parameter region at around Yf- — 10~ 12 and 
8m < 130 MeV is allowed by the observational 7 Li. Here 
5m < 130 MeV corresponds to t? > 10 3 s. The case for 
the destruction of ( 7 Lif~) through the internal conver- 
sion is also similar to that of ( 7 Bef~) (l4l Il5j. 

Further constraints come from the relic density of the 
DM, which can be stated in terms of the stau relic density. 
It is calculated as shown in Fig. 2] for the present values 
of parameters. Applying all the constraints, we are led to 
the allowed interval shown by the thick line in the figure. 



D. Constraint on parameter space of stau 

Finally, we show in Fig.[5]the parameter space in which 
the calculated abundances of the DM and that of the 
light elements are consistent with their values from the 
observations. Here, based on the discussion in previous 
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FIG. 5: Parameter space in which the calculated abundances 
of the DM and that of the light elements are consistent with 
their values from the observations. Parameter region sur- 
rounded by black solid (blue dashed) line shows the allowed 
region from the WMAP observation at the 3a (2a) level [l). 
Red crisscross points show the parameters which is consistent 
with observational abundances for the light elements includ- 
ing 7 Li at 3a level. 



subsection, we took 5m = 100 MeV. Parameter region 
surrounded by black solid (blue dashed) line is allowed 
by the relic abundance of the DM from the WMAP ob- 
servation at the 3cr (2a) level which corresponds to 
0.0913 < n DM h 2 < 0.1285 (0.0963 < fl DM h 2 < 0.1213). 
Red crisscross points show the parameters which are con- 
sistent with the observational abundances for the light el- 
ements including 7 Li, where the observational 7 Li abun- 
dance is yielded by literature [37 1. 

The abundance of the light elements constrains the pa- 
rameter space due to the following reasons. First, the re- 
gion where the stau mass is less than 300 GeV is excluded 
since the relic density becomes too small to destruct 7 Li 
sufficiently. Next, the top-left region of the figure is ex- 
cluded since the lifetime of stau becomes too long and 
hence overproduces 6 Li through the catalyzed fusion 
Indeed, the lifetime of stau gets longer as its mass gets 
heavier due to the small phase space of the final state 
Q. On the other hand, its lifetime gets shorter as the 
left-right mixing angle increases in the present parameter 
space Q. As a result, the final allowed region becomes 
as shown by the red crisscrosses in Fig. [5] 

The black solid curve and the blue dashed curve are 
the constraints from the relic abundance of the DM 
as discussed in subsection IIII Al Note that the relic 
abundance is insensitive to the mass difference dm for 
dm < < m% . Combination of the constraints on the abun- 



dance of the light elements and of the DM strongly re- 
stricts the allowed region and leads to 5m ~ 100 MeV, 
m f = (300 - 400) GeV, and sin(9 T = (0.65-1). In 
Fig. [4] these parameter values correspond to the white 
triangular region below the allowed region (thick line). 
Our model can thus provide a handle to the mixing an- 
gle, which has few experimental signals, once the value 



of mf is determined. 



IV. SUMMARY AND DISCUSSION 

We have studied the evolution of stau number density 
in the MSSM coannihilation scenario, in which the LSP 
and the NLSP are the lightest neutralino and the lighter 
stau, respectively, and have a small mass difference 5m < 
0(1 GeV). In this case, stau can survive until the BBN 
era, and provide additional nucleosynthesis processes. It 
is therefore necessary to see how large the relic density 
of stau is at the BBN era. 

We have shown the Boltzmann equations for the cal- 
culation of the relic number density of stau, and have 
found that the number density of stau continues to evolve 
through the exchange processes Eq. ([TTJ)) . even after the 
relic abundance of the DM is frozen out. Thus, we need 
to calculate the stau relic density by a two-step proce- 
dure. In the first step, we calculate the total abundance 
of stau and neutralino, which corresponds to the relic 
abundance of the DM. The total abundance is controlled 
only by pair annihilation processes of the supersymmet- 
ric particles. In the second step, we calculate the ratio 
of the stau number density to the total number density 
of stau and neutralino, which is governed only by the ex- 
change processes. We have calculated the relic density 
of stau at the BBN era by solving the Boltzmann equa- 
tions numerically. The freeze-out temperature T|( rat ; ) 
is determined by (m T — 5m) /Tf^ at - 1Q ) ~ 25 and the relic 
density of stau are given by Eq. (|23l) . Thus it becomes 
larger as the mass difference between stau and neutralino 
gets smaller. Our method of calculation is generally ap- 
plicable to obtain the density of exotic heavy particles 
that coannihilate with other (quasi)stable particles: we 
calculate the total number density of these particles and 
then calculate the ratio among them by evaluating the 
exchange processes such as Eq. (fTTJ|) . This method of 
calculation can be found versatile in various scenarios in- 
cluding the catalyzed BBN and the exotic cosmological 
structure formation. 

At the BBN era, the long-lived stau form bound states 
with nuclei, and provide exotic nucleosynthesis processes. 
One of them is the internal conversion process, which 
offers a possible solution to the 7 Li problem. Apply- 
ing the calculated relic density of stau, we have calcu- 
lated the primordial abundance of light elements includ- 
ing these exotic processes. We have found the param- 
eter space consistent with both of the calculational re- 
sults and the observations for the relic abundance of the 
DM and the light elements abundance including 7 Li. We 
have shown a prediction for the values of the parameters 
relevant to stau and neutralino, which is shown in Fig. 
[5] Consistency between the theoretical prediction and 
the observational result, both of the DM abundance and 
the light elements abundance requires 5m ~ 100 MeV, 
mf = (300 - 400) GeV, and sin6» r = (0.65 - 1). 
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